/* [wxMaxima batch file version 1] [ DO NOT EDIT BY HAND! ]*/
/* [ Created with wxMaxima version 13.04.2 ] */

/* [wxMaxima: input   start ] */
stat_cond:[1/a^3-1/a^2-1/b^3+1/b^2=0, a+b=R];
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
tot_energy: E[tot]=1/(2*a^2)-1/a+1/(2*b^2)-1/b+1/(a+b);
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
eq_sol:solve(stat_cond, [a,b]);
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
ev(radcan(rhs(tot_energy)),a:eq_sol[1][1],b:eq_sol[1][2]);
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
E1:expand(rhs(%));
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
ev(radcan(rhs(tot_energy)),a:eq_sol[3][1],b:eq_sol[3][2]);
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
E2:expand(rhs(%));
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
wxplot2d([E1,E2], [R,1.8,5], [y,-0.565,-0.465])$
/* [wxMaxima: input   end   ] */

/* Maxima can't load/batch files which end with a comment! */
"Created with wxMaxima"$
